Magnetic phase transitions and magnetoelectric coupling of GdFeOs single crystals 

probed by low-temperature heat transport 
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The low-temperature thermal conductivity (k) of GdFeOa single crystals is found to be strongly 
dependent on magnetic field. The low-field hi{H) curves show two "dips" for H \\ a and only 
one "dip" for H \\ c, with the characteristic fields having good correspondence with the spin-flop 
and the spin-polarization transitions. A remarkable phenomenon is that the subKelvin thermal 
conductivity shows hysteretic behaviors on the history of applying magnetic field, that is, the it{H) 
isotherms measured with field increasing are larger than those with field decreasing. Intriguingly, 
the broad region of magnetic field (~ 0-3 T) showing the irreversibility of heat transport coincides 
with that presenting the ferroelectricity. It is discussed that the irreversible k{H) behaviors are due 
to the phonon scattering by ferroelectric domain walls. This result shows an experimental feature 
that points to the capability of controlling the ferroelectric domain structures by magnetic field in 
multiferroic materials. 

PACS numbers: 66.70.-f, 75.47.-m, 75.50.-y 
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I. INTRODUCTION 

Low-temperature heat transport is an important phys- 
ical property of solids and is useful for probing many 
kinds of elementary excitations, such as phonons, elec- 
trons, magnons and spinous, etc. Thermal conductiv- 
ity (k) is strongly dependent on the statistical laws of 
these excitations and their transport properties, which 
are directly related to the nature of the ground states 
of materials. For example, the temperature dependence 
of very-low-T thermal conductivity can directly show 
the purely phononic transport, the pairing symmetries 
of superconductors, the Fermi-liquid state of metals and 
the nature of spin liquid in quantum magnets.-i^^^ When 
the scattering between different types of quasiparticles 
is significant, the low-T heat transport can also be an 
effective way to detect such couplings; in particular, 
the spin-phonon or magnon-phonon couplings in insu- 
lating magnetic materials are sometimes easily revealed 
by the measurement of the magnetic-field dependence 
of thermal conductivityi^r— In cases of magnetic exci- 
tations either transporting heat or strongly scattering 
phonons, the magnetic phase transitions, including the 
changes of either the ground state or the spin structure, 
can be sensitively probed by the thermal conductivity 
measurements jii"— The purpose of this work is to get 
more information on the magnetic phase transitions and 
spin-phonon coupling of multiferroic materials by study- 
ing the low-T thermal conductivity. 

Multiferroicity is a result of strong coupling between 
magnetic and electric degrees of freedom in insulators 
and has received a lot of research interests because of its 
application usagc i^^'^" It is found that multiferroic ma- 



terials usually present complex H — T phase diagrams 
and multiple magnetic phase transitions, accompanied 
by the drastic changes of electric properties ji^r^^i As a 
result, the low-T heat transport may show peculiar be- 
haviors at these transitions. We choose GdFeOs as a can- 
didate, which has a distorted perovskite structure with 
an orthorhombic unit cell (Pbnm)^ It is known that 
Fe'^^ spins form an AF order along the a axis below 
T^'^ = 661 K with a weak ferromagnetic (WFM) com- 
ponent along the c axis due to the spin canting in the ac 
plane^^i^ which results from the Dzyaloshinskii-Moriya 
interaction^^- The spin structure of Fe'^"'" ions can be ex- 
pressed as GxAyFz in Bertaut's notatiouj^S where Gx, 
Ay and stand for the spin components along the a, b 
and c axes with the NaCl-type, the layered-type and the 
ferromagnetic-type configurations, respectively. On the 
other hand, Gd^^ moments order antiferromagnetically 
along the a axis below T^'' = 2.5 K and show a G^- 
type spin structurc i^'^'^^i'^^ The ferroelectric polarization 
appears below T^"^ and is considered to originate from 
the spin-exchange striction;^^'^'*''^^ more exactly, the in- 
teraction between adjacent Fe^+ and Gd'^^ layers with 
respective G-type AF ordering?- drives Gd^+ ions to dis- 
place along the c axis so as to induce the ferroelectric 
polarization along the c axis. The magnetic-field depen- 
dences of electric polarization (P) were carefully studied 
and discussed to be related to the transitions of mag- 
netic structures. In the case of H \\ a, the simultaneous 
spin-fiop transition of Gd'^"'" moments and reorientation 
of Fe'^"'" spins occur at ^ 0.5 T^^n^ Across this transition, 
the magnetic structure changes from phase I {GxAyF^ 
for Fe'^"'' spins and GxAy for Gd'^^ moments) to phase II 
^Pg3+. FxCyGz] Gd^+: Gz)] thus, the electric polariza- 
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tion originating from the spin exchange striction shows a 
sudden drop. Upon increasing the magnetic field further, 
Gd^+ spins gradually turn to the direction of magnetic 
field and they are completely polarized at 2 T, where 
the magnetic structure changes from phase II to phase 
III (Fe3+: F^CyGz; Gd^+: F^). In the case of || c, 
the magnetic structure changes from phase I to phase IV 
(Fe^+i G^AyF^; Gd^+i F^) at - 2.5 T due to the simple 
spin-polarization transition of Gd^"*" momentsi^ In both 
cases, the electric polarization decreases to zero as long 
as the Gd'^"'" moments are polarized at high fieldjSi 

In this work, we study the low-T thermal conductiv- 
ity of GdFeOa single crystals and find that the magne- 
tothermal conductivity is rather large, indicating a strong 
spin-phonon coupling in this compound. The magnetic- 
field-induced spin-flop (or spin-reorientation) transition 
and spin-polarization transition are detected by k{H) 
isotherms. One peculiarity of the heat transport of 
GdFeOa is that it shows a hysteresis of the subKelvin 
k{H) curves in a broad range of magnetic field. The pos- 
sible origin of this irreversibility is discussed to be due 
to the ferroelectric-domain walls scattering on phonons. 
This result indicates the capability of magnetic field con- 
trolling the ferroelectric domain structures, which is a 
special but understandable phenomenon of multiferroic 
materials. 



II. EXPERIMENTS 

High-quality GdFeOa single crystals are grown by the 
fioating-zone technique in a flow of oxygeni^ The sam- 
ples for thermal conductivity measurements are cut pre- 
cisely along the crystallographic axes with typical dimen- 
sion of 2.5 X 0.6 X 0.15 mm^ after orientated by using 
the x-ray Laue photographs. The a-axis and c-axis ther- 
mal conductivities (/to and Kc) are measured by a conven- 
tional steady-state technique and two different processes: 
(i) using a "one heater, two thermometers" technique in 
a '^He refrigerator and a 14 T magnet at temperature 
regime of 0.3 - 8 K; (ii) using a Chromel-Constantan 
thermocouple in a ''He cryostat for zero-field data above 
4 KfiSii^ In these measurements, the temperature gradi- 
ent is 2-5 % and < 2 % of the sample temperature for 
temperatures below and above 30 K, respectively. The 
specific heat is measured by the relaxation method in 
the temperature range from 0.4 to 30 K using a com- 
mercial physical property measurement system (PPMS, 
Quantum Design). 



III. RESULTS AND DISCUSSION 

Before presenting the heat transport results of GdFeOa 
crystals, we show in Fig. 1 the low-T specific heat data. 
At very low temperatures, a large peak, which is appar- 
ently of magnetic origin, shows up. As can be seen from 
the inset to Fig. 1(a), this low-T peak seems to consist 
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FIG. 1: (Color online) (a) Temperature dependence of the 
specific heat of GdFeOa single crystal below 30 K. The inset 
displays the data below 4 K, which show a sharp peak at 2.2 
K and a shoulder-like feature at ~ 1 K. (b) The specific data 
plotted in C/T vs T^. The dashed line shows the fitting to 
the high-T data by using the formula of phonon specific heat, 
that is, C = /3T^ + ^sT^ + /^rT^ 



of a sharp peak at 2.2 K and a shoulder- like feature (or 
a weak peak) at ~ 1 K. Note that these data essentially 
reproduce those in an earlier report^i However, there is 
one small difference between two sets of data, that is, 
the data in Ref. 'si' showed the sharp peak at 1.47 K and 
the shoulder-like feature at ^ 2 K. It was discussed that 
those two features are originated from the Neel transition 
of Gd"^"'' moments and a Schottky contribution, respec- 
tively. So the data in Fig. 1 indicate that the Neel tem- 
perature of Gd'^"'" moments of our GdFeOa crystal is 2.2 
K, close to the value of a recent report by a susceptibility 
measurement 

It can be seen that the magnetic contributions to the 
specific heat are important only at very low temperatures 
and are likely to be negligible above ~ 12 K, where the 
data in Fig. 1(a) shows a minimum. So one can make an 
estimation of phonon specific heat from the high-T data 
in Fig. 1. It is known that in the temperature range 0.02 
< T/6'd < 0.1 {9d is the Debye temperature), one had 
better use the low-frequency expansion of the Debye func- 
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FIG. 2: (Color online) Temperature dependences of the a- 
axis thermal conductivity of GdFeOa single crystal in zero 
and 14 T magnetic field parallel to the a and c axis. The data 
in 14 T field for two directions are nearly coincident. The 
dashed line indicates the T^'^ dependence. The inset shows 
the temperature dependence of the phonon mean free path / 
divided by the averaged sample width VK in 14 T magnetic 
field. 



tion, C = fiT^+fi^T^ + PjT'^ + where /3, (i^ and (ij are 
temperature-independent coefHcients.— It is found that 
this formula gives a precise fitting to the experimental 
data above 15 K, as shovifn in Fig. 1(b), with the fitting 
parameters /3 = 5.90 x 10""* J/K^'mol, /Ss = 4.52 x 10""^ 
J/K^mol and jSj = -3.96 x 10"!° J/K^mol. Note that at 
very low temperatures, the T^- and T''- terms are neg- 
ligible and the phonon specific heat shows a well-known 
dependence with the coefficient of /3. 
Figure 2 shows the temperature dependences of Ko for 
GdFeOs single crystals in zero and 14 T magnetic field 
parallel to the a and c axis. As far as the zero-field data 
are concerned, they show a typical phonon transport be- 
havior at relatively high temperaturesii In particular, the 
magnitude of the phonon peak at 18 K is as large as 250 
W/Km, which is rather rare in the transition-metal ox- 
ides and indicates weak crystal defects or impurities. At 
low temperatures, however, there are some features show- 
ing a complexity of the phonon transport in this mate- 
rials. First, even at very low temperatures the Ka{T) 
data show a distinct deviation from the law, a sign 
of the phonon boundary-scattering limit, which appar- 
ently indicates the remaining of significant microscopic 
scattering of phonons.^ Actually, the zero-field curve ex- 
hibits a weak kink-like temperature dependence below 2 
K, which is likely due to the magnon-phonon scattering 
since the Gd'^"'" moments order antiferromagnetically be- 
low 2.2 Ki^i2£ The effect of magnetic field on thermal 
conductivity seems to confirm this possibility. When 14 



T magnetic field is applied, the conductivities below 2 K 
become larger and the kink disappears, which clearly in- 
dicates the negative effect of magnons on the heat trans- 
port, considering that at low temperatures the magnons 
can hardly be thermally excited in high field. A T^-^ de- 
pendence, which is very close to the boundary scattering 
limit, indicates that magnetic scattering on phonons is 
almost smeared out in 14 T field. In addition, the effect 
of strong magnetic field on Ka is essentially isotropic for 
H II a and H \\ c. 

It is possible to estimate the mean free path of phonons 
at fow temperatures and to judge whether the phonons 
are free from microscopic scattering at subKelvin tem- 
peratures. The phononic thermal conductivity can be 
expressed by the kinetic formula Kph — ^Cvpl,^ where 
C = f3T^ is phonon specific heat at low temperatures, 
Vp is the average velocity and I is the mean free path of 
phonon. Here /3 = 5.90 x 10^* J/K^mol is obtained from 
the above specific-heat data and Vp = 1930 m/s can be 
estimated from Deybe temperature Qd using the rela- 
tions p = and Qd = ^(^^)^r^ where N 
is the number of molecules per mole and each molecule 
comprises s atoms, V is the volume of crystal and R 
the universal gas constant. So we can calculate I from 
the 14 T k(T) data and compare it with the averaged 
sample width W = ^^JA/tt = 0.361 mm,iiiS, where A 
is the area of cross section. As shown in the inset to 
Fig. 2, the ratio l/W increases with lowering temper- 
ature and becomes close to one at 0.3 K, which means 
that the boundary scattering limit is nearly established 
at such low temperatures. On the other hand, although 
not shown in the figure, the mean free path of phonons in 
zero field is apparently several times smaller than that in 
14 T field, which obviously demonstrates the significance 
of microscopic phonon scattering in zero field. 

The magnon-phonon scattering is further evidenced by 
the magnetic-field dependences of Ka at low tempera- 
tures, as shown in Figs. 3(a) and 3(b). For both H \\ a 
and II c, the HaiH) isotherms show a reduction at low 
fields followed by an enhancement at high fields. At sub- 
Kelvin temperatures, Ka achieves a saturated value in 
high magnetic fields, but the magnitude weakens grad- 
ually with increasing temperature. Note that the high- 
field enhancement of thermal conductivity can be as high 
as ~ 700%, while the strongest suppression at low fields is 
^ 50%, which clearly indicates that the coupling between 
magnons and phonons is pretty strong.— Another re- 
markable feature is that at subKelvin temperatures there 
are two "dips" at low fields (< 4 T) for H \\ a, while 
only one shallow and broad "dip" emerges for _ff || c. 
Moreover, the field-induced suppression of conductivity 
becomes weaker as increasing temperature. At higher 
temperatures above 1 K, another "dip" appears at high 
fields (> 4 T) for both cases and it becomes broader and 
shifts to higher field with increasing temperature, which 
can be clearly seen in Figs. 3(c) and 3(d), suggesting its 
possible origin from phonon scattering by paramagnetic 
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FIG. 3: (Color online) (a,b) Magnetic-field dependences of the 
a-axis thermal conductivity of GdFeOs single crystal. The 
magnetic fields are applied along the a or c axis. All the data 
are measured in the field sweeping-up process after zero-field 
cooling. (c,d) Zoom in of the low-field data of panels (a) and 
(b). 



FIG. 4: (Color online) Low-field Ka{H) isotherms of GdFeOa 
single crystal at subKelvin temperatures with magnetic field 
along the a and c axis. The data shown with solid symbols 
are measured in the ascending field after the sample is cooled 
in zero field, while the open symbols show the data with the 
descending field, as indicated by arrows. 



moments jiiii^ We have also measured the c-axis thermal 
conductivity for both magnetic-field directions (see the 
Appendix), and found that the behavior of k is essentially 
dependent on the direction of magnetic field rather than 
the direction of heat current, which again demonstrates 
the role of magnons in the heat transport in magnetic 
field. 

Figures 3(c) and 3(d) show the details of low-field 
Ka{H) behaviors. As mentioned above, there are clear 
differences in the Ka{H) isotherms between _ff || a and 
H \\ c. When the magnetic field is parallel to the a axis, 
the two "dips" of Ka{H) locate at ~ 0.5 T and - 2.25 T, 
which are weakly temperature dependent. Apparently, 
these two "dips" are directly related to the spin-flop 
transition and the spin polarization of Gd'^^ magnetic 
structurefii as suggested by the low-T electric polariza- 
tion datai^ It was discussed that at 0.5 T the magnetic 
structure changes from phase I to phase II, which is as- 
sociated with a simultaneous occurrence of the Fe'^+ spin 
reorientation and the Gd^^ spin flop both from mainly 
along the a axis to mostly along the c axis. In Bertaut's 
notation, the magnetic structure changes from GxAyFz 
to FxCyGz for the Fe^+ spins and from GxAy to Gz for 
the Gd'^'*' moments, respectively. With increasing mag- 
netic field further, the Gd'^"'" moments gradually rotate to 
the direction of magnetic field and are fully polarized at ^ 
2.25 T, where the magnetic structure changes to phase III 



(represented as F^CyGz for Fe^+ spins and for Gd'^"'" 
moments). Since the magnon excitations become gap- 
less at either the spin-flop fleld or the spin-polarization 
field, the quickly increased number of magnons can scat- 
ter phonons strongly and cause a drastic decrease of ther- 
mal conductivity r^ii^iii^ manifested as those two weakly- 
temperature-dependent "dips" in k[H) curves. When 
the magnetic field is applied along the c axis, there is no 
spin-flop transition of the magnetic structure; instead, 
the Gd"^"*" ions change their spin direction under a simple 
spin polarization process with increasing field. It was 
discussed that the Gd'^"'" ion turns spin direction gradu- 
ally from almost a axis to complete c axis at ~ 2.5 T, and 
the magnetic structure changes from phase I to phase IV 
{GxAyFz for Fe"^"*" spins and Fz for Gd^"*" moments).— 
Consequently, there is only one broad "dip" in Ka{H) for 
H\\c. 

Therefore, all the above results can be well understood 
in the scenario of phonon heat transport, with the signifi- 
cant scattering by magnetic excitations. However, this is 
not the whole story. It is found that the low-T thermal 
conductivity of GdFeOa is dependent on the history of 
applying magnetic field. 

As shown in Fig. 4, for both H \\ a and H \\ c, the 
Ka{H) isotherms measured with field increasing are larger 
than those with field decreasing, forming a clear hystere- 
sis at low temperatures. This phenomenon demonstrates 
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that at very low temperatures, there is some peculiar 
channel of phonon scattering that is related not only to 
the magnetic field but also to the history of applying field. 
As discussed on Figs. 2 and 3, the magnons are effective 
phonon scatterers in zero and low fields. However, they 
cannot simply produce an irreversible behavior of k{H) 
in such a broad field range from almost zero to ~ 3 T. In 
particular, the spin polarization is a naturally continuous 
transition and the hysteretic behavior cannot be expected 
across this transition. On the other hand, other micro- 
scopic phonon scatterers like point defects or dislocations 
in crystal structure are also irrelevant for two reasonsii 
First, the scattering processes from crystal imperfections 
have nothing to do with the external magnetic field, let 
alone the history of applying field; second, these scatter- 
ings are known to be less effective upon lowering temper- 
ature, but the hysteresis becomes more pronounced with 
decreasing temperature. 

Then, what else can be the origin of this kind of scat- 
tering that produces the hysteresis? Considering the de- 
velopment of hysteresis upon lowering temperature, a 
natural origin is related to the magnetic or ferroelec- 
tric domains in this multiferroic material, in which the 
domain walls can play an important role in scattering 
phonons at very low temperatures when the mean free 
path of phonons is long enough to be comparable to the 
interval of domain walls. It is known from the former 
work that the weak-ferromagnetism-related domains in 
GdFeOa can be produced only at very low field (< 0.1 
T)<^ Therefore, the ferroelectric domain walls are likely 
the main source of phonon scattering that is responsible 
for the k{H) irreversibility in the broad field range well 
above 0.1 T. In this regard, there are several experimen- 
tal results supporting the origin of k{H) hysteresis from 
the ferroelectric domains. First, the electric polarization 
vs magnetic field P{H) also showed the irreversible be- 
haviors in the absence of electric field for both H \\ a 
and H \\ c.~ Second, the field range that the hysteretic 
behavior of n{H) appears is almost the same as that ex- 
hibiting the hysteretic P{H) curvesi^ Third, the rela- 
tive magnitude of polarization for field increasing and 
field decreasing indicated that there are less ferroelec- 
tric domains when the field is increasing .-^^ It is there- 
fore in good agreement with the present observation that 
the thermal conductivity is larger in the field-increasing 
branch, which can be due to the weaker phonon scat- 
tering by domain walls. Furthermore, with increasing 
temperature the size of hysteresis diminishes and disap- 
pears completely above ~ 1 K, also in good agreement 
with the fact that the phonon scattering by boundaries 
are unimportant at high temperatures where the mean 
free path of phonons becomes much shorter than the av- 
eraged distance of domain walls. 

The magnetic-field-induced irreversibility of k can also 
be manifested in the k{T) curves, as shown in Fig. 5, 
measured in some characteristic fields by warming up 
the sample after zero-field cooling or field cooling to the 
lowest temperature. It is found that below 1 K the zero- 
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FIG. 5: (Color online) Temperature dependences of the a- 
axis thermal conductivity of GdFeOa single crystal in 2 T 
magnetic field along the a axis (a) and in 2.5 T field along 
the c axis (b). The data are taken with slowly warming up 
the sample from the lowest temperature after field cooling or 
zero-field cooling. 



field cooled (ZFC) conductivities are larger than the field 
cooled (FC) ones, which coincides with the hysteresis of 
n{H) curves. This difference indicates that there are less 
ferroelectric domains if the sample is cooled in zero field. 
From these data, one can get a rough estimation of the 
averaged distance between the ferroelectric domain walls, 
which usually has the same order of magnitude to that 
of the mean free path of phonons as long as the domain 
walls are effective phonon scatterers. Following the above 
calculation on the 14 T k{T) data, the mean free path of 
phonons is obtained to be 20-30 /xm from the ZFC data 
in Fig. 5. This is a rather reasonable value of domain 
size for the ferroelectric single crystals. 

Note that the hysteresis of n{H) of a magnetic material 
itself is not strange. The irreversibility has been known 
to often appear at the field-induced first-order magnetic 
transition, such as the spin-Peierls to AF order— and the 
liquid-gas-like transition in spin-ice compounds^ etc. 
However, the hysteretic k{H) of GdFeOs is somewhat 
different, because it appears in a rather broad field range 
where the magnetization does not show irreversibility)^ 
In contrast, another well-known large hysteresis of low- 
T n{H) is the one observed in the high- Tc cuprate 
BizSrsCaCusOg (BSCCO)^i28 in which the hysteresis 
is well understood by the vortex-pinning effect that is 
accompanied with the irreversible macroscopic magneti- 
zation. From the above discussions, one can see that the 
irreversible k{H) of GdFeOs is most likely caused by the 
ferroelectric domain walls scattering on phonons, that is 
dependent on the history of applying the magnetic field. 
It should be pointed out that the phonon scattering by 
the ferroelectric domain walls had been studied by mod- 
ifying the domain structure through applying the elec- 
tric field directly i^^i^" It was proved in KH2PO4 that the 
low-T thermal conductivity of the single-domain state is 
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much larger than that of the multi-domain statej^ Sim- 
ilar result was obtained in SrTiOs, in which the low-T 
thermal conductivity under high electric field is clearly 
larger than that in zero fieldf^ However, these former 
experiments should be analyzed very carefully because 
applying an electric field using the contacts on the sam- 
ple surface could bring some uncertainties in the ther- 
mal conductivity measurements. In this regard, chang- 
ing the domain structures by applying magnetic field is 
free from such problem and can give more reliable data. 
It is intriguing that in multiferroic GdFeOa the ferroelec- 
tric domain structures are able to be manipulated by the 
magnetic field, which is a key to make the peculiar low-T 
k{H) hysteresis to be observable. 



IV. SUMMARY 

The low-temperature heat transport of GdFeOs single 
crystal is found to be strongly dependent on the mag- 
netic field and is helpful for studying the magnetic transi- 
tions and magnetoelectric coupling in this material. The 
magnons play an important role in the low-T thermal 
conductivity by scattering phonons instead of acting as 
heat carriers. As a result, both the spin-flop or reorienta- 
tion of the magnetic structure and spin polarization cause 
significant phonon scattering at the transition fields. The 
most remarkable result is that the low-T (< 1 K) ther- 
mal conductivity shows an irreversible behavior on the 
history of applying magnetic field. This phenomenon is 
rather peculiar in the sense that it appears in a broad 
range of magnetic field (~ 0-3 T) where the magneti- 
zation is known to be almost reversible. Furthermore, 
the hysteresis of k{H) or the irreversibility of ZFC and 
FC k(T) become larger upon lowering temperatures. All 
these data suggest that the ferroelectric domain struc- 
tures, which are surprisingly sensitive to the history of 
applying magnetic field, are playing the key role in pro- 
ducing the irreversibility of thermal conductivity through 
the phonon scattering by the domain walls. The present 
results show an interesting case that the ferroelectric do- 
mains of a multiferroic material can be manipulated by 
magnetic field, and this kind of magnetoelectric coupling 
can be observed by the heat transport measurement at 
very low temperatures. 
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FIG. 6: (Color online) Temperature dependences of the c-axis 
thermal conductivity of GdFeOs single crystal in zero and 
14 T magnetic field along the a or c axis. The inset shows 
the temperature dependence of the phonon mean free path I 
divided by the averaged sample width W (= 0.344 mm) in 14 
T magnetic field. 



Appendix: The c-axis thermal conductivity of 
GdFeOa 



Figure 6 shows the temperature dependences of Kc for 
GdFeOs single crystals in zero and 14 T magnetic field 
parallel to the a and c axis. Apparently, these data are 
essentially the same as those of Ka{T), including all the 
main features like the magnitude of phonon peak, the 
weak-kink-like temperature dependence below ~ 2 K, the 
significant recovery of conductivity in 14 T field, etc. In 
addition, the temperature dependences of 14 T thermal 
conductivities are also similar to those of Ka{T) in 14 
T. The calculated mean free path of phonons along the 
c axis is found to have the same order of magnitude to 
that from Ka, as shown in the inset to Fig. 6. One can 
easily conclude that the phonon heat transport is nearly 
isotropic in GdFeOs. 

Figures 7(a) and 7(b) show the magnetic-field depen- 
dence of the c-axis thermal conductivity of GdFeOs single 
crystal. The overall behaviors of Kc{H) are qualitatively 
the same as those of Ka{H), that is, the thermal conduc- 
tivity are suppressed at low fields and strongly enhanced 
at high fields. Moreover, at low fields there are also two 
"dips" for H || a and one "dip" for || c, respectively. 
One may note that the "dip" fields of Kc for both H \\ a 
and II c are somewhat different from those of Ka when 
the field is applied along the same direction. However, 
such discrepancy is likely due to the demagnetization ef- 
fect. Both samples have size about 2.5 x 0.6 x 0.15 
mm'^. For such kind of long-shaped samples, the demag- 
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FIG. 7: (Color online) (a,b) Magnetic-field dependences of 
the c-axis thermal conductivity of GdFeOa single crystal. The 
magnetic fields are applied along the a or c axis and all the 
data are measured in the field sweeping-up process after zero- 
field cooling, (c-f) Low-field Kc{H) isotherms of GdFeOs 
single crystal at subKelvin temperatures with magnetic field 
along the a and c axis. The solid symbols show data measured 
with ascending field after the sample is cooled in zero field, 
while the open symbols show the data with descending field, 
as indicated by arrows. 



netization factor n is negligible (« 0) when the applied 
field is along the longest dimension, while it is not negli- 
gible (taking a value between and 1) when the field is 
along the shortest dimension. For the a-axis sample, its 
a axis is along the longest dimension and the c axis along 
the shortest dimension; while the c-axis sample has the 
c axis and a axis along the longest and the shortest di- 
mensions, respectively. Therefore, for the experimental 
configurations of (kq, H \\ a) and (kc, H \\ c), the in- 
ner magnetic field (Hi) is equal to the external field (H) 
since n w 0; while for the configurations of («;„, H |j c) 
and (kc, H \\ a), the inner field is different to the exter- 



E 




0.01 



T(K) 



T(K) 



FIG. 8: (Color online) Temperature dependences of the c- 
axis thermal conductivity of GdFeOa single crystal in 2.5 T 
magnetic field along the a axis (a) and in 2 T field along 
the c axis (b). The data are taken with slowly warming up 
the sample from the lowest temperature after field cooling or 
zero-field cooling. 



nal one and can be expressed as Hi — H/{l + nx), where 
X is the magnetic susceptibility. Note that although the 
factor (1 -f- nx) is not available for our samples, it is pos- 
sible to make a comparison between Ka{H) and Kc{H) 
data with taking into account the demagnetization ef- 
fect. The comparison of the "dip" fields between Ka{H) 
and Kc{H) for both H \\ a and H \\ c indicated that the 
factor (1 -I- nx) is about 1.5. For example, in the case 
of {Ka, H II a) the "dip" fields of 0.5 T and 2.25 T are 
intrinsic values of transition fields; while the "dip" fields 
of 0.75 T and 3.25 T in the case of (kc, H \\ a) re-scaled 
by the factor of 1.5 give values of 0.5 T and 2.2 T, respec- 
tively, which both match the "dip" fields of Ka{H) very 
well. Under such kind of consideration, it is a bit strange 
that the "dip" fields of (kq, H \\ c) and {kc, H \\ c) do not 
differ much from each other. Therefore, the quantitative 
analysis of these data are called for. 

The Kc{H) isotherms for both H \\ a and H \\ c also 
present the hysteresis on the magnetic field sweeping up 
and down, as shown in Figs. 7(c-f). It can be seen 
that the field region and temperature regime for showing 
the hysteresis, the difference between the field up and 
down data are all very similar to those of Ka{H) data. 
Apparently, these results have the same origin as those 

of Ka{H). 

Figure 8 shows the magnetic-field-induced irreversibil- 
ity of Kc(T), measured in some characteristic fields by 
warming up the sample after zero-field cooling or field 
cooling to the lowest temperature. They also show simi- 
lar behaviors to those of the ZFC and FC Kq(T) data. 

In summary, the low-T heat transport of GdFeOa is 
nearly isotropic. There are strong coupling between 
the phonons and magnons, which determines the main 
features of the magnetic-field dependence of thermal 
conductivity. The characteristic transitions of k(H) 
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isotherms are strongly dependent on the direction of ex- 
ternal field and are caused by the field-induced transi- 
tions of magnetic structure. An irreversibility of low-T 
thermal conductivity on the magnetic field is observed 
and is likely due to the phonon scattering by the ferro- 



electric domain walls. The peculiarity is that this phe- 
nomenon points to the manipulation of ferroelectric prop- 
erties by magnetic field, which is a characteristic of mul- 
tiferroic materials. 
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